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simplificada de un índice de sequía, empleando una 
distribución de probabilidad mezclada. Tecnología y Ciencias 
del Agua, 7(5), 135-149. 


En la actualidad, muchos son los procedimientos para 
caracterizar el fenómeno de las sequías, la mayoría de ellos 
se basa en el Índice Estandarizado de Sequía (SPI), el cual 
se calcula a partir de registros históricos de precipitación 
para un cierto periodo de estudio. Esta serie de datos 
se estandariza empleando diversas distribuciones de 
probabilidad, lo que permite el análisis de los impactos de 
las sequías a distintas escalas de tiempo. Lo complejo de 
este procedimiento es el uso de diversas distribuciones de 
probabilidad y que ninguno de sus parámetros tiene un 
significado físico. Este artículo presenta una metodología 
para el cálculo simplificado del índice SPIL con base 
en la hipótesis de que una serie de eventos ocurre de 
forma aleatoria en el tiempo. Se emplea un proceso de 
Poisson para explicar el número de días de lluvia y una 
distribución exponencial para caracterizar las láminas 
medias de precipitación. Buscando una explicación física 
de los parámetros de dicha distribución de probabilidad 
mezclada (Poisson-exponencial) se correlacionan de modo 
directo con el índice SPI, obteniendo un índice modificado 
que simplemente se calcula contando el número de días 
de lluvia del periodo, representando de manera sencilla y 
confiable la condición de sequía. La metodología se aplica 
en ocho estados del norte de la república mexicana. 


Palabras clave: sequía, distribución de Poisson, régimen de 
precipitación, índice estandarizado de precipitación. 


Abstract 


Gutiérrez-López, A., Vargas-Baecheler, J., Reséndiz-Torres, V., 
é Cruz-Paz, 1. (September-October, 2016). Simplified Design 
of a Drought Index, Using a Mixed Probability Function. Water 
Technology and Sciences (in Spanish), 7(5), 135-149. 


Nowadays, there are many procedures to characterize the 
phenomenon of drought; most of them are based on the Standardized 
Precipitation Index (SPD), which is calculated based on historical 
records of rainfall for a certain period of study. This set of data is 
standardized using various probability distributions, which enables 
the analysis of the impacts of droughts at different time scales. 
The complexity of this process is the use of various probability 
distributions and that none ofits parameters have a physical meaning. 
This paper presents a simplified methodology for calculating the SPI, 
based on the hypothesis that a series of events occur randomly in 
time. A Poisson process is used to explain the number of raining 
days and an exponential distribution to characterize the mean 
rainfall high. Looking for a physical explanation of the parameters 
of the distribution of mixed probability (Poisson-exponential), they 
correlate directly with the SPI index, obtaining a modified index is 
calculated simply by counting the number of raining days in the 
period; representing a simple and trustworthy way the condition of 
drought. The methodology is applied in eight northern states of the 
Mexican Republic. 


Keywords: Drought, Poisson distribution, rainfall regimen, 
standardized precipitation index. 
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Introducción 


En la hidrología moderna es común observar 
series de eventos hidrometeorológicos que ocu- 
rren de manera aleatoria en el tiempo (Govin- 
daraju, 2013). Los procesos de Poisson propor- 
cionan modelos probabilísticos apropiados para 
explicar el número de eventos que ocurren o se 
producen a partir de cierto fenómeno (López- 
Segovia, Villaseñor-Alva, € Vaquera-Huerta, 
2002). En hidrología es frecuente relacionar la 
distribución de Poisson con la ocurrencia de 
los eventos de precipitación, para explicar, por 
ejemplo, el intervalo de ocurrencia de tormentas 
como variables discretas, en donde se evalúa la 
ocurrencia o no de un evento de lluvia; intensi- 
dad y duración se presentan como dos variables 
aleatorias independientes y discretas (Tapsoba, 
1997; Bacchi, Becciu, éz Kottegoda, 1994; Le 
Barbé, Lebel, éz Tapsoba, 2002; Gutiérrez-López, 
Descroix, éz Lebel, 2002a). Incluso estudios tan 
precisos, como cuantificar el número de gotas, 
se basan en esta distribución. Calder (1986) 
propone un modelo estocástico de modelación 
de precipitación, el cual relaciona, mediante la 
distribución de Poisson, el número promedio 
de gotas de lluvia retenidas en una superficie 
dada con el número promedio de gotas de lluvia 
caídas. Lo anterior fundamenta la afirmación de 
que el proceso de Poisson es una herramienta 
adecuada para representar el comportamiento 
de fenómenos hidrometeorológicos (Ribstein, 
1983; Wong, Lambert, Leonard, €: Metcalfe, 2010; 
Sol áková De Michele, €: Vezzoli, 2014; Sohrabi, 
Ryu, Abatzoglou, €: Tracy, 2015). Cuando se 
utiliza el modelo mezclado Poisson-exponencial 
es posible evaluar también el riesgo hidrológico 
asociado con inundaciones y sequías de forma 
simultánea (Gutiérrez-López, Lebel, € Descroix, 
2002b) de forma similar a un proceso jerarquiza- 
do de Neyman-Scott (Rodríguez-Iturbe, 1986). 
En esta distribución mezclada se demuestra que 
los parámetros de las dos distribuciones están 
relacionados a través de un factor de reducción 
de varianza, el cual es función de la estructura 
de correlación espacial de la precipitación (Cam- 
malleri, Micale, £z Vogt, 2016). Suponiendo que 
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la lámina media por evento de precipitación si- 
gue una distribución exponencial, el parámetro 
de escala ajustado a partir de la serie de lluvia 
diaria adquiere el mismo valor del parámetro 
de escala de la distribución de valores extremos 
tipo Gumbel, ajustando los valores máximos 
de lluvia mensual (Sivapalan € Blóschl, 1998; 
Gutiérrez-López, Lebel, € Mejía, 2005). Tra- 
bajos similares reemplazando la distribución 
exponencial por la distribución gamma han 
sido desarrollados en cuencas con registros 
de precipitación horaria en el centro del Reino 
Unido y Estados Unidos (Eagleson, 1978; Onof 
é£x Wheater, 1994; Wilks, 1998). Asimismo, la 
distribución de Poisson suele emplearse para 
modelar periodos de sequía (Sharma, 1996; 
Cameron, Beven, éz Tawn, 2000). Por lo general, 
los periodos de sequía y de humedad durante 
una temporada de lluvias tienden a conservarse 
y pueden ser modelados utilizando un proceso 
de Markov de orden uno. Una forma tradicional 
para el análisis de estos fenómenos mínimos 
consiste en mezclar más de dos distribuciones. 
Así, una función de probabilidad de Poisson 
delimita la ocurrencia de periodos de sequía; la 
distribución geométrica se utiliza para obtener 
las probabilidades asociadas con el número de 
periodos, y la distribución Weibull modela la 
precipitación total dentro de un periodo húme- 
do (Mirakbari, Ganji, € Fallah, 2010). El análisis 
completo, que requiere la estimación conjunta 
de los parámetros de estas distribuciones, arroja 
resultados relacionados con la probabilidad de 
que en un día cualquiera se presente sequía; la 
probabilidad de que un día seco sea consecutivo 
a Otro día seco, o la probabilidad de que un día 
húmedo sea consecutivo a otro día húmedo. Si 
bien este tipo de resultados describen de forma 
adecuada el intervalo de tiempo entre periodos 
de sequía, es deseable correlacionar los pará- 
metros de una distribución de probabilidad, 
por ejemplo, con algún índice de aridez. Éste es 
precisamente el tercer objetivo de este trabajo. 
Para el caso de caudales en un río, un esquema 
similar de detección del déficit o exceso de 
corrientes en un intervalo de tiempo puede 
ser modelado por una distribución de Poisson 
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(Coles, Pericchi, €: Sisson, 2003; Hughes, 2003). 
Debido a que en la mayoría de los países los 
ingenieros civiles evalúan las crecientes de 
diseño utilizando la extrapolación de valores 
extremos para el diseño de estructuras hidráu- 
licas, el pronóstico hidrológico se realiza a 
partir de una interpretación de las frecuencias 
de probabilidades de ocurrencia de los eventos 
máximos (Pandey € Nguyen, 1999). En vir- 
tud de que el nombre de valores extremos es 
asignado a estas distribuciones debido a que 
éstas pueden ser obtenidas como distribuciones 
límite de los valores grandes o pequeños entre 
variables aleatorias independientes, cada una 
teniendo una distribución continua (Pandey éz 
Nguyen, 1999; Escalante $z Reyes, 2002), puede 
entonces hacerse la hipótesis de que el proceso 
compuesto Poisson-exponencial representa fe- 
nómenos extremos de sequía o de inundación 
debido a que contiene una distribución de ex- 
tremos. Esto permitiría asociar la distribución 
mezclada Poisson-exponencial con una ley de 
valores extremos; es decir, eventos definidos por 
un proceso de Poisson, que dan lugar a lluvias 
puntuales distribuidas de forma exponencial. 


Proceso compuesto Poisson-exponencial 
(Ley Fuga) 


Esta distribución mezclada de probabilidad tie- 
ne su origen en el estudio de series cronológicas 
de eventos separados por duraciones aleatorias. 
Las dos hipótesis de base son la estacionariedad, 
que significa que las propiedades estadísticas de 
la serie cronológica no cambian con el tiempo, y 
la independencia, que se refiere a que cada uno 
de dichos eventos ocurre de manera indepen- 
diente de los demás. El nombre de “Ley Fuga” 
(Loi de Fuites) proviene del estudio que realizó 
M. Morlat (Babusiaux, 1969) en la distribución 
de los gastos que se fugaban de una red de dis- 
tribución de gas en Francia. Consiste en utilizar 
simultáneamente una distribución de Poisson, 
para representar la probabilidad de ocurrencia 
de un evento y una distribución exponencial pa- 
ra expresar la magnitud media de dicho evento. 


Si se considera una variable aleatoria y > 0, 
que sigue una distribución exponencial simple, 
con parámetro de escala f, su función de densi- 
dad puede expresarse como: 


sum= he (1) 


Sea £ un número entero aleatorio que sigue 
una Ley de Poisson, de parámetro de forma 
A > 0, el valor de su probabilidad será: 


en 


ñ (2) 


p,(U)= 


Así, una variable aleatoria X = E 
Y, puede representarse como la suma de 2 > 0 
variables aleatorias independientes que siguen 
un proceso compuesto Poisson-exponencial con 
parámetro de forma h, que denota la tasa de lle- 
gadas de un evento por unidad de tiempo, y un 
parámetro de escala f$. Respecto a la duración, es 
equivalente decir que un número de eventos n(t) 
seguirá una Ley de Possion para una duración 
t, O bien que el intervalo T entre dos eventos 
consecutivos sigue una ley exponencial sobre 
el mismo tiempo t. 


£0= 0 (8) 


donde: 
A: duración entre eventos. 
T: intervalo entre dos eventos. 


Cabe recordar que la distribución exponen- 
cial es un caso particular de la llamada familia 
de distribuciones Gamma (Saporta, 1990). Sea 
t,=T,+T, +...+ T, la suma de r variables que 
siguen una distribución exponencial (ecuación 
(3)). Es sabido que la función de densidad de la 
suma de r variables aleatorias independientes 
es el producto de convolución de las densidades 
de cada una de ellas (Tapsoba, 1997), lo cual per- 
mite considerar T como una variable que define 
las probabilidades de tener 1, 2, o más eventos 
desde cero hasta T. 

Así, la función de densidad de probabilidad 
de la llamada Ley Fuga (proceso compuesto 
Poisson-exponencial) para la variable reducida 
U= X/f estará dada por la expresión: 
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fee” Led (4) 
[O=f=e* 


Donde: 


A es el parámetro de forma A > 0 (número de 
eventos, días con lluvia dentro de un periodo 
de tiempo 7). 

PB  esel parámetro de escala f¿ > 0 (lámina media 
de lluvia por evento). 

I, es la función de Bessel de primer orden. 


Es importante mencionar que este proceso 
compuesto Poisson-exponencial tiene una dis- 
continuidad para el valor cero (excepto cuando 
A = 1); es decir, la ordenada al origen tiene un 
valor de he”. La figura 1 muestra las funciones 
de distribución y probabilidad de dicho proceso. 


Estimación de parámetros 


Considerando las hipótesis de esta distribu- 
ción de probabilidad y el significado de los 
parámetros, se infiere que un procedimiento 
simple es el número de eventos nulos n, que 
existan en la muestra; esto es, el número de días 
secos, en la serie histórica de datos. Si n,%0, el 
procedimiento de estimación requiere, para una 
correcta solución, tomar en cuenta el tamaño n 
de la muestra (Tapsoba, 1997). El logaritmo de 
la función de verosimilitud, de la cual se busca 
obtener el máximo, es la siguiente: 


jaa E Y Log,x, 
2 Pp i=1 B 2 i21 


EAS (5) 


La anulación de las derivadas parciales de 
L respecto a A y f conducen a un sistema de 
ecuaciones: 


sE Y A =0 (6) 
24 Ei J24 
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n == Mo ng X. 1n-Ng L(z,) YA 
z + i_ 2 =( (7) 
24 2 $ 216328 
en donde: 
A La! (8) 
B 


[,(z) es la primera derivada de 1,(z). Resol- 
viendo el sistema de ecuaciones miembro a 
miembro de las ecuaciones (6) y (7) se obtiene: 


21 

Il 

> 
pa) 


(9) 


Las propiedades de la función de Bessel 
permiten escribir la derivada de [, como: 


L2)=1,6)-21,(2) 


[, es la función de Bessel de orden cero. 
Empleando esta expresión en la ecuación (5) y 


sustituyendo X en la expresión (9) > = 1% se 
A 


] ye 
obtiene: 


ny 


10) Y Vin mO 


cuya solución es un valor único y positivo de A. 
La función g(A) es monótona y decreciente, su 


solución se encuentra en el intervalo [107 mom) 


en donde A corresponde al valor del paráme- 
tro estimado por el método de momentos. El 
valor de f se obtiene de la expresión (8). Le Bar- 
bé y Lebel (1997) proponen un procedimiento 
alternativo para la estimación de los parámetros 
del proceso compuesto Poisson-exponencial. 
Este método propone emplear el número 
promedio de días secos n, dentro del periodo 
T del total de n días. Tal procedimiento tiene 
una singular importancia, pues la información 
de los días secos no se emplea de manera directa 
en la estimación de parámetros de los métodos 
de momentos. De esta forma se tiene que los 
parámetros estará dados por: 


h--tog| 2. b-> (10) 
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Figura 1. Funciones de densidad y acumulada para diferentes valores de A. 


U=X/P 


Aplicación de la Ley Fuga 


Para ejemplificar la aplicación de esta dis- 
tribución en la modelación y análisis de la 
variabilidad espacial de la precipitación se 
seleccionaron ocho estados del norte de México, 
empleando series de lluvia diaria en registros 


promedio de 50 años de extensión dentro de 
un periodo que comprende los años de 1960 a 
2013. El número de estaciones es variable, sin 
embargo se considera bien representado para 
cada uno de los estados seleccionados: Aguas- 
calientes (49 estaciones); Estado de México (38); 
Guanajuato (21); Jalisco (85); Michoacán (20); 
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Querétaro (31); San Luis Potosí (28), y Zacatecas 
(29). Para la validación del modelo se seleccionó 
la Región Hidrológica 10, situada al noroeste de 
la república mexicana. Las lluvias en esta zona 
se presentan en una época bien definida, pues 
el periodo más abundante se registra entre los 
meses de junio y octubre (Conagua, 2014). Sin 
embargo, es también una región sensible y ex- 
puesta a los cambios climáticos y a las sequías. 
Se han registrado largos periodos de sequía y 
también inundaciones provocadas principal- 
mente por fenómenos meteorológicos extremos. 


Validación del modelo 


Se utilizaron las series históricas de las 93 
estaciones climatológicas dentro de la región 
hidrológica 10 y se estimaron los parámetros 
del proceso compuesto Poisson-exponencial, 
empleando los tres procedimientos ya descritos. 
Con base en el estadístico x? se seleccionaron 
los mejores ajustes y el método de momentos 
demostró, en la mayoría de los casos, dar el 
mejor ajuste a las muestras de datos (Gutiérrez- 
López et al., 2002b). Como se mencionó antes, 
el parámetro f ofrece una ventaja adicional, ya 
que con éste puede deducirse la distribución de 
la lluvia considerando diferentes intervalos de 
tiempo. Así, los parámetros A* y $* deducidos en 
forma simplificada pueden estimarse de forma 
directa multiplicando por el intervalo de tiem- 
po deseado. Es decir, los parámetros tienen un 
significado físico específico y pueden ser usados 
para diferentes valores de T. De esta forma, se 
tiene que el parámetro f$ permanece constante 
a lo largo del tiempo y el parámetro de forma 
A puede calcularse para cualquier duración T, 
esto es: 


B=B m=TA, (11) 


Se aceptará que existe no estacionaridad en 
el caso de que el valor de A tenga una variación 
mayor a 20% entre el valor máximo y mínimo 
dentro del periodo considerado. Esto quiere de- 
cir que se considerará como un proceso estacio- 
nario mientras se mantenga constante durante 
80% del periodo; en caso contrario, el parámetro 
de escala debe estimarse como (Gutiérrez-López 
et al., 2002b): 


ficas (12) 


Cuando se dispone de datos de lluvia diaria 
es recomendable elegir T igual a un día (Lebel 
éz Le Barbé, 1997; Lebel, Taupin, £ D'Amato, 
1997). El cuadro 1 muestra la validación del 
modelo para la estación "Tamazula para el 
periodo 1947-2005. Los parámetros estimados 
del análisis de la lluvia diaria por el método de 
momentos fueron 4 = 0.654 y K =12.214. De ellos 
se pueden deducir los parámetros, por ejem- 
plo, para un intervalo de 10 días, empleando 
la ecuación (4). Esto es: A* = 0.654 x 10 = 6.54 
y PB* = 12.214. Mientras que los valores reales 
de los parámetros, estimados de la serie de 
datos original, para el intervalo de 10 días son 
h= 7.071; P= 11.465. De la misma forma se pue- 
den estimar los parámetros para la cualquier 
otra duración T deseada. 

Una ventaja de utilizar este proceso 
compuesto Poisson-exponencial radica en la 
posibilidad de asociar sus parámetros con 
un significado físico, los cuales pueden ser 
relacionados con conceptos de la gestión del 
riesgo hidrológico, ya que el parámetro $ puede 
asociarse con una medida de la probabilidad 
de lluvia extrema y "3 puede relacionarse con la 
probabilidad de ocurrencia de lluvia (sequía), 
tal y como se desarrollará a continuación. 


Cuadro 1. Validación de los valores de los parámetros de escala y forma del proceso compuesto Poisson-exponencial. 


Estimado de la serie de datos original 


Estimado con la ecuación (4) 


t=1 día A=0.654 y $=12.214 


t=10 días 2=7.071 y B = 11.465 WM =6.54 y P* = 12.214 
t=30 días 2. = 20.083 y B =12.33 AM =20.274 y P*= 12.214 
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El parámetro A de la Ley Fuga como 
índice de sequía 


El índice de precipitación estandarizado SPI 
(Standardized Precipitation Index) es uno de 
los índices utilizados en el mundo para eva- 
luar y monitorear las sequías (Palmer, 1965; 
McKee, Doesken, €: Kleist, 1993; Edwards é 
McKee, 1997). En general, este índice se calcula 
para intervalos de tiempo entre 3 (valor de SPI 
mínimo) y 48 meses (valor de SPI máximo), y 
permite identificar y clasificar posibles zonas 
susceptibles a presentar algún tipo de daño por 
sequía. La metodología para el cálculo del SPI se 
fundamenta en el ajuste de una serie de registros 
históricos de precipitación total mensual a una 
función de distribución Gamma incompleta, y 
en la transformación de los datos resultantes 
a la función de distribución normal estándar. 
Cabe recordar que la deducción del proceso 
compuesto Poisson-exponencial empleando una 
distribución exponencial es un caso particular 
de la llamada familia de distribuciones Gamma. 
De esta forma, el SPI representa el número de 
desviaciones estándar que cada registro de 
precipitación se desvía del promedio histórico. 
Los registros de precipitación superiores al 
promedio histórico del mes correspondiente 
tendrán valores positivos del SPI (presencia de 
humedad), mientras que los registros de precipi- 
tación inferiores al promedio histórico arrojarán 
valores del SPI negativos (déficit de humedad), 


tal como se muestra en el cuadro 2. Si bien se han 
propuesto otros índices de sequía simplificados 
(Pita, 2003; Karamouz, Rasouli, e Nazif, 2009), 
a la fecha ningún trabajo ha incluido de modo 
explícito el número de días lluviosos (intuiti- 
vamente incluido en la precipitación mensual) 
para estimar un índice de sequía. La aplicación 
del proceso compuesto Poisson-exponencial 
ofrece una ventaja, pues establece una relación 
del parámetro de posición de dicha distribución 
de probabilidad con el índice SPI (Kim, Valdés, 
éz Yoo, 2006; Liu et al., 2013). Basados en el hecho 
de que un evento de sequía puede caracterizarse 
definiendo su régimen, es decir, su magnitud 
y frecuencia, parece apropiado introducir la 
variable correspondiente al número de días con 
lluvia. 

De esta forma, se propone relacionar el índi- 
ce tradicional de sequía SPI con el parámetro A 
del proceso compuesto Poisson-exponencial. A 
partir de la estimación de parámetros propuesta 
por Babusiaux (1969), empleando la ecuación 
(10) y poniendo todo en función del índice 
SPI, se obtiene una relación SPI-A, que permite 
estimar un índice modificado SPI*. Esta relación 
es de la forma: 


y2 
SPI*(A)= 0-2 a E - Log, > 
n S, Mo 


+Log, (13) 
n 


0 


Cuadro 2. Valores convencionales de los índices SPI y PDSI(NDMC, 2003; Hayes, 2000). 


SPI Condición PDSI 

+2.00 o más Excepcionalmente húmedo +5.00 a +6.00 
+1.50 a +2.00 Extremadamente húmedo +4.00 a +5.00 
+1.20 a +1.50 Muy húmedo +3.00 a +4.00 
+0.80 a +1.20 Moderadamente húmedo +2.00 a +3.00 
+0.50 a +0.80 Humedad incipiente +1.00 a +2.00 
+0.50 a -0.50 Humedad cercana a la normal +1.00 a -1.00 
-0.5 a -0.8 Sequía incipiente -1.00 a -2.00 
-0.8 a -1.2 Sequía moderada -2.00 a -3.00 
-1.2a-15 Sequía severa -3.00 a -4.00 
-1.50 a -2.00 Sequía extrema -4.00 a -5.00 
-2.00 o menos Sequía excepcional -5.00 a -6.00 
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Donde: 


A  esel número de eventos (días con lluvia) 
dentro de un periodo de tiempo T. 

SPI* es el índice de sequía estimado a partir del 
número de días con lluvia. 


La figura 2 muestra la relación que existe 
entre el índice SPI y el valor del parámetro A 
de la distribución de probabilidad Poisson- 
exponencial, representado con la letra L, para 
el mes de enero en el estado de Querétaro. De 
la misma forma, se obtuvieron las correlaciones 
para el resto de los estados analizados. En todos 
los casos se observa que existe una correlación 
promedio de 0.90, lo que garantiza una estima- 
ción confiable de la relación que existe entre el 
número de días con lluvia (A) y el índice SPI. 
Las ecuaciones que predominan son los ajustes 
a rectas con valores agrupados en el valor de 
A < 0.40 y después valores dispersos para valores 
superiores. A continuación se evalúan todos los 
meses del año para los estados en estudio. Los 


resultados de este procedimiento se muestran 
en la ilustración 3. En este caso y por razones de 
espacio, sólo se muestran algunos de los ajustes 
más representativos en los estados, en los meses 
de mayor sequía histórica. 


Discusión 


Una vez que se verifica la correlación del índice 
de sequía SPI y el valor del parámetro A de la 
distribución de probabilidad Poisson-exponen- 
cial para todos los meses más representativos de 
los estados analizados, se procede a calcular el 
valor de SPI* como un estimado del valor real 
de SPI. Es importante mencionar que en algunos 
de los estados de nuestro país, la condición del 
SPI medio es muy cercana a la condición del SPI 
mínimo, como puede observarse en las figuras 
3b, 3d y 3e. 

En estas condiciones ocurre en los meses 
de invierno (diciembre, enero y febrero) lo 
que caracteriza el régimen de precipitaciones 
típico de la zona central de nuestro país, con 
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Figura 2. Correlación del índice de sequía SPI y el valor del parámetro A del proceso compuesto Poisson-exponencial, para el 


mes de enero en el estado de Querétaro. 
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Figura 3. Correlación del índice de sequía SPI y el valor del parámetro A del proceso compuesto Poisson-exponencial, para 


algunos de los meses más representativos de los estados analizados 
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escasa precipitación en dicho periodo. Para un 
mes históricamente seco, como abril, la figura 
3f muestra las condiciones mínima y máxima 
también muy próximas. Sin embargo, para los 
meses considerados dentro de la temporada ex- 
tendida de lluvia (mayo y noviembre), inicio y 
final de la temporada, respectivamente (figuras 
3a y 3c), los patrones en el ajuste de los valores 
de SPI y el valor del parámetro A se observan 
bien definidos, y separadas las condiciones 
mínima y media. A continuación se estiman 
los valores del índice SPI* para cada uno de 
los estados seleccionados, en cada una de sus 
estaciones climatológicas. Estos resultados se 
muestran en los cuadros 3 y 4, en forma com- 


parativa con el valor del SPI tradicionalmente 
calculado según McKee et al. (1993), y Edwards 
y McKee (1997). Se puede observar cómo en 
todos los casos, excepto uno, la condición de 
sequía es la misma. Cabe mencionar que si bien 
el valor calculado no es el mismo, una vez que 
se compara la condición de sequía resulta la 
misma (cuadro 2). Los resultados muestran 
que el cálculo del parámetro SPI* preserva la 
condición de sequía de un evento, tal que éste 
se inicia cuando se observa a través del tiempo, 
una tendencia marcada de valores negativos 
continuos del SPI. De esta manera, el evento de 
sequía termina cuando el valor del SPl alcanza 
valores positivos. Sin embargo, como puede in- 


Cuadro 3. Condición de sequía para el estado de Querétaro, según la formulación tradicional del índice SPI y la comparación 
con el índice SPI* basado en el número de días de lluvia del mes de enero. 


SPI 
(McKee et al., 1993; (Poi o iS 
a E Edwards de McKee, 1997) oisson-exponencial, 2016) 
sPI A A A o E 
media) Clasificación de la sequía SPI* = 3.364L - 0.5303 Clasificación de la sequía 

22001 [|Amealco, Amealco 0.195 Humedad cercana a la normal 0.121 Humedad cercana a la normal 

22002 |Ayutla, Arroyo Seco 0.198 Humedad cercana a la normal 0.121 Humedad cercana a la normal 

22003 [Cadereyta, Cadereyta 0.274 Humedad cercana a la normal 0.121 Humedad cercana a la normal 

22004 |El Batán, Corregidora 0.292 Humedad cercana a la normal 0.229 Humedad cercana a la normal 

22006 [El Pueblito, Corregidora 0.310 Humedad cercana a la normal 0.338 Humedad cercana a la normal 
= 22008 |Jalpan, Jalpan 0.320 Humedad cercana a la normal 0.338 Humedad cercana a la normal 
E 22011 [Pedro Escobedo 0.342 Humedad cercana a la normal 0.338 Humedad cercana a la normal 
E 22012 |Peñamiller Peñamiller 0.356 Humedad cercana a la normal 0.338 Humedad cercana a la normal 
3 22015 [San Juan del Río (SMN) 0.359 Humedad cercana a la normal 0.446 Humedad cercana a la normal 
= 22017 |Tolimán, Tolimán 0.362 Humedad cercana a la normal 0.446 Humedad cercana a la normal 
E 22018 [Villa Corregidora 0.384 Humedad cercana a la normal 0.446 Humedad cercana a la normal 
o 22029 |Huimilpan, Huimilpan 0.484 Humedad cercana a la normal 0.555 Humedad cercana a la normal 
E 22031 [Presa Paso de Tablas, Teq. 0.502 Humedad cercana a la normal 0.555 Humedad cercana a la normal 
E 22032 |Presa Constitución, San Juan 0.506 Humedad cercana a la normal 0.555 Humedad cercana a la normal 
5 22033 [San Pablo, Amealco 0.514 Humedad cercana a la normal 0.555 Humedad cercana a la normal 
E 22036 |Arroyo Seco, Arroyo Seco 0.550 Humedad cercana a la normal 0.555 Humedad cercana a la normal 
= 22042 |La Llave, San Juan del Río 0.555 Humedad cercana a la normal 0.555 Humedad cercana a la normal 
? 22045 |Juriquilla, Querétaro 0.571 Humedad cercana a la normal 0.663 Humedad cercana a la normal 
= 22046 [Nogales, Colón 0.572 Humedad cercana a la normal 0.663 Humedad cercana a la normal 
S 22047 [Presa El Capulín, Amealco 0.637 Humedad cercana a la normal 0.663 Humedad cercana a la normal 
S 22049 |El Zamorano, Colón 0.707 Humedad cercana a la normal 0.663 Humedad cercana a la normal 
E 22051 [|El Comedero, Peñamiller 0.873 Humedad cercana a la normal 0.880 Humedad cercana a la normal 
S 22056 |Higuerillas, Cadereyta 1.391 Muy húmedo 1.097 Moderadamente húmedo 
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Cuadro 4. Condición de sequía para el estado de Querétaro, según la formulación tradicional del índice SPI y la comparación 
con el índice SPI* basado en el número de días de lluvia del mes de mayo. 


ao 1993; a 
a E Edwards 8: McKee, 1997) (Poisson-Exponencial, 2016) 
sPI Aa 4 ri 4 
míñimo) Clasificación de la sequía |SPI*=2,1631 Ln(L) + 1.1158| Clasificación de la sequía 
22001 [Amealco, Amealco -2.070 Sequía excepcional -2.103 Sequía excepcional 
22002 |Ayutla, Arroyo Seco -1.720 Sequía extrema -1.332 Sequía severa 
22003 [Cadereyta, Cadereyta -1.650 Sequía extrema -1.332 Sequía severa 
22004 |El Batán, Corregidora -1.500 Sequía extrema -1.125 Sequía moderada 
22006 |El Pueblito, Corregidora -1.090 Sequía moderada -1.125 Sequía moderada 
22008 |Jalpan, Jalpan -1.050 Sequía moderada -1.125 Sequía moderada 
22011 [Pedro Escobedo -1.020 Sequía moderada -1.125 Sequía moderada 
22012 |Peñamiller, Peñamiller -0.970 Sequía moderada -1.125 Sequía moderada 
22015 [San Juan del Río (SMN) -0.840 Sequía moderada -1.125 Sequía moderada 
22017 |Tolimán, Tolimán -0.840 Sequía moderada -1.054 Sequía moderada 
22018 [Villa Corregidora -0.820 Sequía moderada -0.937 Sequía moderada 
22029 | Huimilpan, Huimilpan -0.810 Sequía moderada -0.764 Sequía incipiente 
22031 |Presa Paso de Tablas, Teq. -0.810 Sequía moderada -0.764 Sequía incipiente 
22032 |Presa Constitución, San Juan | -0.780 Sequía incipiente -0.764 Sequía incipiente 
22033 |San Pablo, Amealco -0.750 Sequía incipiente -0.764 Sequía incipiente 
22036 | Arroyo Seco, Arroyo Seco -0.750 Sequía incipiente -0.764 Sequía incipiente 
22042 |La Llave, San Juan del Río -0.750 Sequía incipiente -0.764 Sequía incipiente 
22045 |Juriquilla, Querétaro -0.730 Sequía incipiente -0.764 Sequía incipiente 
22046 |Nogales, Colón -0.670 Sequía incipiente -0.764 Sequía incipiente 
22047 |Presa El Capulín, Amealco -0.660 Sequía incipiente -0.604 Sequía incipiente 
22049 |El Zamorano, Colón -0.610 Sequía incipiente -0.604 Sequía incipiente 
22051 |El Comedero, Peñamiller -0.580 Sequía incipiente -0.604 Sequía incipiente 
22056 |Higuerillas, Cadereyta -0.480 | Humedad cercana a la normal -0.454 Humedad cercana a la normal 


tuirse, para una estimación exhaustiva de una 
condición de sequía o aridez se deben estimar 
otros índices empleando distintas variables, 
como precipitación, evapotranspiración, 
caudales en cauces naturales, niveles de lagos 
y embalses, niveles de acuíferos, etc. (Maity, 


Sharma, Nagesh-Kumar, €£ Chanda, 2013). 
Conclusiones 


La lluvia en el norte de la república mexi- 
cana es sin duda un fenómeno complejo de 
estudiar. Sin embargo, los resultados de la 
aplicación del proceso compuesto Poisson- 
exponencial en esta zona del país han permi- 


tido un conocimiento simple del fenómeno, 
a partir de sólo contar el número de días de 
lluvia. La Ley Fuga, como se le conoce, mues- 
tra ventajas. En un principio facilitó conocer la 
distribución espacial y temporal de la lluvia, a 
través de la obtención de los dos parámetros 
de la distribución, tal como lo demostraron 
Chulsang, Daeha, Tae-Woong y Kyu-Nam 
(2008). La sencillez de dichos parámetros 
permitió encontrar relaciones simples con su 
entorno geográfico y físico. 

Se demuestra que el modelo propuesto es 
significativamente consistente con la descrip- 
ción de los cambios asociados con los fenómenos 
extremos, pues permite identificar en una forma 
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confiable la condición de sequía de una zona a 
través de la evolución del parámetro A. Para la 
mejor aplicación del procedimiento propuesto 
se presenta el Anexo A, en donde se muestran 
el total de ecuaciones para los ocho estados 
analizados a lo largo del año. Las ecuaciones son 
muy sencillas y representan la variabilidad del 
fenómeno de sequía de 1960 a 2013. Es suficiente 
tener acceso a datos diarios de precipitación y 
realizar un análisis de frecuencia con la ecuación 
(4), estimando sus parámetros con las ecuacio- 
nes (5) a (10). A continuación se debe verificar 
la propiedad aditiva de los parámetros, según 
se propone en las ecuaciones (11) y (12). Por 
último, se cuenta el número de días con lluvia 
para el mes analizado; por ejemplo, seis días de 
lluvia promedio para el mes de enero en Ameal- 
co, Querétaro, esto es, A = 6/31 = 0.1936, que 
sustituyendo este valor en la ecuación de SPI”* 
media presentada en el Anexo A, se obtiene SPI* 
= 0.121, condición de sequía equivalente a “hu- 
medad cercana a la normal”; misma condición 


obtenida de la aplicación del cálculo tradicional 
del índice SPI. Es así como se considera que la 
metodología propuesta y el empleo del proceso 
compuesto Poisson-exponencial se adaptan a las 
necesidades de estimación del riesgo hidrológi- 
co y permiten profundizar en el conocimiento 
de la climatología regional del norte de la repú- 
blica mexicana. 

Respecto a la comparación entre el índice SPI 
y el número de días con lluvia, se puede obser- 
var que el parámetro A por sí solo representa 
de forma aceptable un riesgo de sequía y que 
permite obtener a través de una relación sencilla 
el índice tradicional de sequía SPLI. El proceso 
compuesto Poisson-exponencial permite esti- 
mar de manera rápida y confiable un índice de 
sequía que en lo general es difícil estimar y sobre 
todo de interpretar. Cabe mencionar que estos 
índices de aridez suelen apreciarse también a 
través de diversas características fisiográficas o 
climatológicas de la región, como el potencial 
de evaporación (Arora, 2002). 


Anexo A. Ecuaciones para el cálculo de la condición de sequía (SPD), empleando el 


número de días de lluvia A. 


Aguascalientes Estado de México Guanajuato Jalisco 
Mín. |SPI=2.1198 L - 1.067 SPI = 4.7697 L - 1.7198 SPI = 1.1772In(L)+1.0521 SPI = 5.5132 L + 0.0094 
SS Enero Med. |SPI=0.9012 L - 0.0738 SPI = 3.6879 L221% SPI = 1.4797L - 0.2155 SPI = 14.028 L - 0.9193 
Ne Máx. | SPI=5.0045 L+ 0.9168 SPI = 4.4137 L05%7 SPI = 1.1794In(L) + 3.742 SPI = 23.337 L + 1.3916 
0) 
E Mín. —|SPI=0.7793In(L) +1.1632 | SPI = 4.4369 L - 1.5511 SPI= 1.00511n(L) +1.5377 SPI = 3.3928 L - 0.5588 
> Febrero Med. |SPI=0.392In(L) + 0.9858 | SPI=2.4675 L'”*2 SPI=0.46161n(L) + 1.111 SAL = 7/7 ILLES 
E Máx.  |SPI= 6.7212 L + 0.7142 SM= 1 A ls SPI= 1.4956In(L) +4.7404 SPI= 3.3779 L + 1.5388 
o 
a Mín. —|SPI=3.6147 L- 0.1466 SPI = 1.1938 In(L)+0.9374 | SPI=9.0688L - 1.9433 SPI = 0.6636 In(LW-1.7166 
Y Marzo Med. |SPI=2.5172 L + 0.2527 SPI = 2.5389 2215 SPI = 15.265L2?%9 SPI = 0.5052 In(L)+1.6211 
3 Máx. |SPI=7.0078 L + 1.1569 SPI = 4.4326 L03125 SPI = 15.769L17 SPI = 0.8193 In(L)+4.0796 
ja] 
2 Mín. |SPI=2.8415 L - 0.6699 SPI = 4.7189 L - 3.3246 SPI= 4.0933L - 1.6687 SPI = 3.9557 L - 0.6343 
8 Abril Med. |SPI=1.541 L02 SPI =0.002e77%4L SPI=1.3071L - 0.1653 SPI = 2.7022 L - 0.055 
3 Máx. |SPI=3.726 L + 1.2808 Ser 0:93752: 925 SPI=7.8405L + 0.1171 SPI= 1.9066 e! 1298L 
E Mín. —|SPI=1.1697In(L) + 0.578 | SPI=3.4139 In(L) -0.0745 | SPI=4.1397L - 2.7999 SPI = 4.19 L - 1.9891 
5 Mayo Med. |SPI= 0.8373 L - 0.1355 SPI = 0.4564 L - 0.2397 SPI=0.6292L - 0.1684 SPI = 1.8807 L?1%7+ 
E Máx. |SPI=6.6314 L + 0.2172 SPI = 2.2894 In(L) + 2.93 SPI = 4.3293L + 0.3558 SPI = 1.4133 e *ert 
E 
= Mín. —|SPI=2.574In(L)- 0.03 SPI = 3.8822 L - 4.6468 SPI = 8.0486L - 7.1462 SPI = 3.6454 In(L) -0.1636 
E Junio Med. |SPI= 0.4708 L - 0.2181 SPI = 0.4265 L - 0.2984 SPI=0.8645L - 0.5467 SPI= 0.3736 L - 0.2054 
E Máx.  |SPI= 3.8053 L - 0.1456 SPI = 0.4736€e!081 1 SPI = 5.433L - 1.9128 SPI = 0.5287 e20056L 
3, Mín. —|SPI=3.756 L - 4.4262 SPI = 5.3854 L - 6.0263 SPI = 4.2992L - 4.8967 SPI = 6.7904 In(L) -0.4224 
= Julio Med. |SPI=0.4633 L - 0.32 SPI = 0.4932 L - 0.3623 SPI=0.461L - 0.3249 SPI = 0.4645 L - 0.3521 
.S Máx. | SPI=4.4924 L'% SPI = 2.8201 2036 SPI = 3.6083L12 SPI =0.1669e*'75L 
$ Mín. —|SPI=2.7258 In(L) -0.2203 | SPI = 4.4066 L - 5.0641 SPI = 4.3371L - 4.6674 SPI = 3.8731 In(L) -0.6868 
y) Agosto Med. |SPI=0.3582 L - 0.1731 SPI = 0.4884 L - 0.3521 SPI = 0.4874L - 0.3086 SPI= 0.2748 L - 0.1795 
2 Máx.  |SPI=3.7586 L128 Ser 023982025 SPI =2.0734In(L) +2.7367 SMI= MAL eto 
50 
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Mín. SPI = 2.6318 In(L)+0.1261 | SPI = 3.5855 In(L) -0.7107 | SPI=4.1565L - 4.5179 SPI = 2.758 In(L) + 0.9552 
Septiembre | Med. |SPI=0.5114 L - 0.2193 SPI = 0.4904 L - 0.3476 SPI = 0.3862L72! SPI = 1.1181 L - 0.4077 
Máx. —|SPI=3.5437 L185% SPI = 3.4947 In(L)+2.8464 | SPI=2.9811L!'"”> SPI = 4.0283 L + 0.5607 
Mín. SPI = 1.4517 In(L) + 0.506 | SPI = 1.562 In(L) - 0.3399 SPI = 1.4285In(L) +0.0708 SPI = 1.1693 In(L) -0.3281 
Octubre Med. |SPI=0.9981 L?20 SPI = 0.3521 L - 0.121 SPI = 0.2338In(L) +0.2811 SPI = 0.1354 In(L)+-0.1809 
Máx.  |SPI=3.7774 L+ 0.6306 SPI = 1.5364 In(L)+2.7965  |SPI=1.1913In(L) +3.0004 SPI=2.8101 L*72 
Mín. SPI = 3.217 L - 0.9645 SPI = 1.4547 In(L)+0.7965 | SPI =0.8591In(L) +0.9083 SPI = 0.5719 In(L)+0.3922 
Noviembre | Med. |SPI=1.38 L - 0.0104 SPI = 1.003 L - 0.2136 SPI = 1.3032L - 0.0627 SPI = 0.24 In(L) + 0.5619 
Máx. | SPI= 4.4486 L0188 SPI = 2.1622 In(L)+4.9147 | SPI=1.1702In(L) +4.0784 SPI = 0.608 In(L) + 3.0718 
Mín. SPI = 4.0484 L - 1.1414 SPI = 6.1976 L - 1.9903 SPI=0.5906In(L) +0.6796 SPI = 0.4932In(L) +0.3827 
Diciembre | Med. |SPI=1.883 L - 0.1188 SPI= 6.5009 L21748 SPI=0.261n(L) + 0.6988 SPI=0.1783In(L) +0.5155 
Máx. |SPI=4.1526 L0*%2 SPI = 1.1834 e2é8L SPI = 1.2157In(L) +4.2605 SPI = 4.7815L + 0.9552 


Anexo A. Ecuaciones para el cálculo de la condición de sequía (SPI), empleando el 
número de días de lluvia A (continuación). 


Michoacán Querétaro San Luis Potosí Zacatecas 

Mín. SPI = 18.513L - 1.118 SPI = 5.3415 L - 1.4758 SPI = 2.4119 In(L) + 1.5291 SPI = 2.4463 L - 1.0427 
Enero Med. |SPI=8.6146L - 0.115 SPI = 3.364 L - 0.5303 SPI = 1.1419 L - 0.2781 SPI = 0.0921 e*929L 

Máx. SPI = 15.218L+ 1.7685 SPI = 4.4982 L + 0.9289 SPI = 1.1103 In(L) + 3.5595 SPI = 3.1158 L + 1.5664 

Mín. SPI = 14.639L - 0.6033 SPI = 1.2802 In(L) + 2.0643 | SPI = 1.8446 In(L) + 1.5617 SPI = 0.6648 In(L) + 1.022 
Febrero Med. |SPI=8.0164L + 0.1014 SDI = 0.6274 In(L) + 1.428 | SPI = 0.3972 In(L) + 0.7029 SPI = 0.3599 In(L) + 0.9581 

Máx. SPI = 15.897L + 1.6192 SPI=5.8657 L99% SPI = 0.9402 In(L) + 3.6397 SPI = 1.1444 In(L) + 4.0017 

Mín. SPI = 0.5979In(L) +2.4724 | SPI= 0.9164 In(L) + 1.596 | SPI=2.1855 In(L) + 2.1661 SPI = 0.842 In(L) + 2.0056 
Marzo Med. |SPI=12.293L + 0.1617 SPI = 0.5067 In(L) + 1.2898 | SPI =2.187 L - 0.4009 SPI = 3.6414 L + 0.0453 

Máx. SPI = 23.199L + 1.7119 SPI = 5.3723 L + 1.0837 SPI = 3.5904 L + 1.2149 SPI = 1.3176 In(L) + 4.9408 

Mín. SPI= 6.176L + 0.0648 SPI = 5.2691 L - 2.0731 SPI = 6.0426 L - 2.8676 SPI = 3.2475 L - 0.5198 
Abril Med. |SPI=11.265L - 0.6988 SPI = 4.5624 L?8192 SPI = 1.2722 L - 0.2726 SPI = 1.5818 L + 0.1081 

Máx. SPI = 11.981L + 1.9738 SPI = 5.2516 L + 0.3947 SPI = 3.1819 L + 1.2224 SPI = 3.4021 L + 1.6052 

Mín. SPI = 5.3881L - 1.7338 SPI = 2.1631 In(L) + 1.1158 | SPI =1.5886 In(L) + 0.007 SPI = 1.4083 In(L) + 1.1657 
Mayo Med. |SPI=1.1792L - 0.0489 SPI= 1.1099 L - 0.311 SPI = 0.4205 L - 0.123 SPI = 1.22 L - 0.1845 

Máx. SPI =7.0514L + 1.1123 SPI = 3.1182 L + 0.6043 SPI = 2.405 L + 0.9533 SPI = 1.2642! 916L 


Mín. SPI=0.6164L - 4.217 SPI = 4.1222 L - 3.4613 SPI = 4.557 L - 3.9842 SPI = 2.2934 In(L) + 0.0825 
Junio Med. SPI = 0.0577L - 0.1919 SPI = 0.7549 L - 0.3397 SPI = 0.4759 L - 0.194 SPI = 0.5374 L - 0.2291 
Máx. SPI = 0.4764L - 0.053 SPI = 2.8141 L + 0.3001 SPI = 3.3646 L + 0.3275 SPI = 3.0751 L + 0.3572 
Mín. SPI = 5.6131L - 4.7579 SPI = 2.346 In(L) + 0.428 SPI = 3.4218 L - 3.5257 SPI = 2.6702 L - 3.062 
Julio Med. SPI = 0.4956L - 0.2288 SPI = 0.8609 L - 0.3636 SPI = 0.3977 L - 0.1827 SPI = 0.3136 L - 0.1519 
Máx. SPI = 2.5988L + 0.5955 SPI = 3.2038 L0%%48 SPI = 2.5079 L + 0.691 SPI = 3.2977 L + 0,33 
Mín. SPI = 4.4194L - 3.8375 SPI= 2.1166 In(L) + 0.4269 | SPI=3.0105L - 3.007 SPI = 2.1946 In(L) - 0.1036 
Agosto Med. SPI = 0.4546L - 0.1947 SPI = 0.7152 L - 0.2544 SPI = 0.3189L - 0.1069 SPI = 0.4546 L - 0.2063 
Máx. SPI = 4.2408L' SPI = 3.1514 L + 0.2581 SPI = 1.0944e!2827L SPI = 3.1946 L0882 
da Mín. SPI = 2.9536ln(L) +0.8106 | SPI = 3.6749 L - 3.066 SPI = 2.8075 L - 3.127 SPI = 1.7118 In(L) - 0.1787 
eptiembre 
P Med. SPI = 0.6012L - 0.225 SPI = 0.8371 L?9063 SPI = 0.3212 L - 0.1294 SPI = 0.4399 L - 0.1749 
Máx. SPI = 5.0473L1-14% SPI = 3.4711 L - 0.0141 SPI = 2.5281 L + 0.692 SPI = 1.3817 In(L) + 2.9076 
Mín. SPI = 1.9056ln(L) +1.5694 | SPI = 1.487 In(L) + 0.7302 | SPI= 3.7443 L - 2.8359 SPI = 3.0948 L - 2.0169 
Octubre Med. SPI = 0.9724L - 0.1336 SPI= 0.3526 In(L) + 0.5312 | SPI= 0.7301 L?737 SPI = 0.7654 L17 
Máx. SPI = 1.2109In(L) + 3.904 | SPI=2.9386 L + 0.7578 SPI=5.1117 L=0658 SPI = 4.6949 L + 0.4538 
Mín. SPI =0.77811n(L) +1.5378 | SPI=1.0301 In(L) + 1.3444 | SPI = 1.4483 In(L) + 0.6895 SPI = 2.6008 L - 0.9177 
Noviembre | Med. SPI = 0.2507In(L) +0.8839 | SPI = 0.3818 In(L) + 0.9088 | SPI = 1.8983 L?*5 SPI = 1.2978 L - 0.0566 
Máx. SPI =0.6232In(L) +3.9355 | SPI = 3.0197 L%218 SPI = 0.9799 In(L) + 3.6353 SPI = 3.4942 L + 1.3534 
Mín. SPI=0.695In(L) + 1.8087 | SPI=0.9637 In(L) + 1.6571 | SPI= 1.5359 In(L) + 1.1605 SPI = 2.6381 L - 1.0126 
Diciembre | Med. SPI = 5.3223L - 0.0089 SPI= 0.5204 In(L) + 1.3119 | SPI= 0.3778 In(L) + 0.6587 Sie = (ula Y Mt 
Máx. SPI = 18.057L + 1.4411 SPI = 0.952 In(L) + 3.807 SPI = 1.0866 In(L) + 4.022 SPI = 4.088 L0306 
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